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Abstract 

We find that the fermion determinant is not generally positive in a class 
of lattice actions recently constructed by Cohen et al. hep-lat/0307012 ; these 
are actions that contain an exact lattice supersymmetry and have as their 
target (continuum) theory (4,4) 2-dimensional super- Yang-Mills. We discuss 
the implications of this finding for lattice simulations and give some preliminary 
results for the phase of the determinant in the phase-quenched ensemble. 
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1 Introduction 



A major motivation for efforts to latticize supersymmetric models is that some non- 
perturbative aspects of supersymmetric field theories are not accessible by the usual 
techniques, most often relying on holomorphy. Indeed, phenomenological applica- 
tions of softly broken A/" = 1 4d supersymmetric gauge theories such as the minimal 
supersymmetric standard model are generally regarded as effective field theories valid 
for the TeV regime, derived from a more fundamental supersymmetric theory valid 
at higher energy scales. To give rise to spontaneous supersymmetry breaking, the 
more fundamental theory must contain other sectors that are responsible for this 
effect. Often the additional sectors include a strongly- interacting super- Yang-Mills 
theory. The most general effective theory consistent with symmetry constraints typ- 
ically involves nonholomorphic quantities, determined by the Kdhler potential. It is 
of interest to understand the nonperturbative corrections that this potential may re- 
ceive, particularly if the corresponding dynamics are expected to play an important 
role in determining qualitative features of the effective theory; in addition to breaking 
supersymmetry, nonperturbative effects are generally considered to be important to 
moduli stabilization (vacuum selection). 

Thus, one hope for lattice supersymmetry is that it would lead to simulations that 
would provide further data on nonperturbative aspects of supersymmetric field the- 
ories, especially those that include super- Yang-Mills. However, the lattice regulator 
generically breaks supersymmetry. This is a natural consequence of the fact that the 
supersymmetry algebra is embedded into the super-Poincare algebra, which involves 
both rotations and translations. Only a discrete subgroup of the Poincare group 
survives, so it is not surprising that the regulator is not supersymmetric. Thus the 
target theory is obtained by fine-tuning the bare parameters of the lattice action — 
supplemented by counterterms; a considerable amount of work has been done in this 
direction. Early efforts in this direction can be found in Pl E] ; recent reviews with ex- 
tensive references are given in O E] ■ In the case of pure 4d A/" = 1 super- Yang-Mills, 
a gauge invariant formulation with chiral lattice fermions suffices to guarantee the 
correct continuum limit; for example, simulations using domain wall fermions have 
been performed 

It is worth noting, however, that the discrete rotations and translations that are 
preserved in a typical (isotropic, hypercubic) latticization guarantee that in the con- 
tinuum limit the only relevant and marginal operators that can appear are those that 
preserve the full Poincare invariance of the target theory. It is of interest to explore 
whether or not it is possible to realize something analogous for the super-Poincare 
group. Some amount of exact lattice supersymmetry — i.e., a fermionic symmetry that 
relates lattice bosons and lattice fermions and closes with the discrete subgroup of 
the Poincare group that is preserved by the lattice — would presumably allow for the 
target theory to be obtained in a more controlled fashion. Remarkably, in some recent 



1 



examples possessing an exact lattice supersymmetry, it has been found that certain 
super- Yang-Mills theories may be obtained from a lattice theory without the need for 
fine-tuning; additionally, other types of super- Yang-Mills theories may be obtained 
with less fine-tuning than would occur in a naive discretization 13 E] ■ It is worth 
noting that other sorts of models with exact lattice supersymmetry have been dis- 
cussed in the literature by a few groups: supersymmetric quantum mechanics P ITU], 
the 2d Wess-Zumino model ^T] pure super- Yang-Mills using overlap fermions 
[T^ , and direct constructions in the spirit of the Ginsparg- Wilson relation^ have all 
been considered. In some cases, the constructions have been related to topological 
field theory [17\ . 




1.1 Deconstruction, orbifolds and lattice supersymmetry 

In this letter we will be interested in the super- Yang-Mills constructions that lead 
to a Euchdean lattice whose target theory contains 8 supercharges jSj- This is a 
generalization of the 4 supercharge constructions of [7] . The method of building all 
such models is based on deconstruction of extra dimensions [IHllini- The correspond- 
ing interpretation in terms of the world-volume theory of D-branes has led to the 
latticizations of 2d, 3d and 4d supersymmetric gauge theories. While a complete 
latticization of spacetime is studied here, partial latticizations also yield interesting 
results; for example, the chiral anomalies [201 or instanton solutions [21] of 4d = 2 
super- Yang-Mills can be equivalently described by deconstructing one dimension to 
obtain a 3d = 4 product group (quiver) super- Yang-Mills theory that is an effective 
latticization of one dimension. 

The Euclidean spacetime lattice constructions are all arrived at by orbifold pro- 
jections of super- Yang- Mills matrix models.^ I.e., in each case we quotient a matrix 
model by some discrete symmetry group of the theory. Degrees of freedom that are 
not invariant with respect to the combined action of the orbifold generators are pro- 
jected out.^ Following jHlEj, we will refer to the "nonorbifolded" matrix model as 
the mother theory and the "orbifolded" matrix model as the daughter theory. The 
(effective) lattice theory is obtained by studying the daughter theory expanded about 
a nontrivial minimum of its scalar potential; i.e., a point in its moduli space. 

^Luscher first suggested that the perfect action approach might be used to identify an analogue 
of the Ginsparg- Wilson relation for lattice supersymmetry |13| . This idea was worked out in nonin- 
teracting examples in |14l I15| . (An approach very similar to |15j has been applied in |lf>j , yielding 
slightly different expressions.) 

^These matrix models are obtained as Od reductions of 4d, 6d and lOd J\f = 1 Euclidean super- 
Yang-Mills. 

•^For a detailed discussion, we refer the reader to [HI [7] ■ 
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1.2 The fermion determinant 



In lattice theories containing fermions, it is well-known that it is of great practical 
importance that the fermion determinant (or more generally, the Pfaffian), obtained 
by integrating over the fermion degrees of freedom in the partition function, be pos- 
itive semi-definite. For let v represent the lattice bosons. Then having integrated 
out the lattice fermions, one obtains an equivalent effective lattice action {Sb is the 
bosonic part of the action before integrating out fermions): 



A positive semi-definite det M{v) yields a real effective action, thus avoiding the 
inherent problems of a complex action with respect to estimating correlation functions 
by Monte Carlo simulation. 

However, in [22j we found that the fermion determinant was complex for the 
constructions of which have (2,2) 2d U{k) super- Yang-Mills as their target theory.^ 
Here we study whether or not a similar problem exists in the constructions of [Hj, 
which have (4,4) 2d U{k) super- Yang-Mills as their target theory. 

1.3 Summary of results 

For the reader's convenience, we now summarize the content of our work: 

• In this letter we show that the lattice theory with (4,4) 2d super- Yang-Mills as 
its target obtained from orbifolded supersymmetric matrix models, possesses 
a problematic fermion determinant. 

• Due to a ever-present fermion zeromode, detM(0) = 0, i.e., for all boson con- 
figurations. The zero eigenvalue can be factored out in a controlled way in 
order to exhibit the determinant for the other fermion modes. In the daughter 
theory with (4,4) 2d SU{2) super- Yang-Mills as its target, we carry out this 
factorization (numerically). 

• Once the zeromode fermion has been factored out, we find that the remaining 
product of eigenvalues is generically nonzero with arbitrary complex phase. 

• We discuss the implications of our results for lattice simulations of the latticized 
(4,4) 2d super- Yang-Mills theories. We emphasize that the complex action may 
be turned into a virtue, in that it provides an interesting system in which to 
study complex action simulation techniques — relevant to lattice QCD at finite 
temperature and/or baryon density. 

^In 122 we also showed that this was not in conflict with the well-known positivity of the fermion 
determinant in the mother theory | 23l I24| . 
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• We conclude by presenting some preliminary Monte Carlo results for the dis- 
tribution of the phase of the fermion determinant when sampled in the phase- 
quenched distribution; that is, when the determinant in (^Q) is replaced by its 
absolute value. We discuss why this may be relevant to the continuum limit. 
Unfortunately, for the small lattice that we consider (2x2), we find that the 
phase distribution is essentially uniform in it range (— 7r,7r]. 

In the remainder of this letter we will discuss various details related to these results. 



2 Construction 



2.1 Mother theory 

The action of the mother theory is that of a 6d ^ Od reduction of U{kN'^) JV = 1 
sup er- Yang- Mills : 

^ = Tr i[v^, Vn] Vn]) + ^eij Tr (^f CS„[t;^, ^,]) (2) 

where e = iaa, Vm = v^T"', = *f {i = 1, 2) with = (1, T") a Hermitian basis 
for the generators of U^kN"^). Each of the two is a 4-component fermion. The 
4x4 matrices are components in the construction of the 6d (Euclidean) Clifford 
algebra, and C is a charge conjugation matrix. For further details we refer the reader 
to 0. 



2.2 Orbifold to daughter theory 

In addition to the U{kN'^) gauge invariance, the mother theory possesses an SO{Q)e 
Euclidean invariance group and an SU{2)b. chiral R-symmetry group. From the 
invariance group SO{Q)e ® SU{2)r one isolates a f/(l)i ® U{1)2 subgroup, with 
corresponding generators ri and r2. The bosons and the fermions of the mother 
theory are written in a basis with well-defined ri,r2 charges. In the construction of 
[H] this basis is:^ 

*f = (A, 6, 6, 6), = (x, ^1,^2, ^3), (3) 



/o 



-X 

-y 



X 



—z 

y 



y 

z 


-X 



-y 

X 

0/ 



(4) 



'Our lattice boson notation is related to that of ^ by (a;, y, z) = (zi, Z2, z^). 
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The ri,r2 charges are given in Table 1 of [HI; we will not need them here. All that is 
important is that the bosons and fermions (denoted collectively by $) of the mother 
theory are subjected to a projection with respect to a Zj^ ® Zp^ subgroup of f/(l)i ® 
[7(1)2, together with a nontrivial embeddeding into the gauge group U{kN'^). That 
is, we keep only fields that satisfy 

$ = e^-"''--!^ Ca^C^^ (5) 

where Ca are generators of a Zj^ ® Z^ subgroup of U^kN"^). This breaks the gauge 
group down to U{k)'^^ = <S>mi m2=i ^(^)mi,m25 corresponding the U{k) gauge invari- 
ance of an x lattice theory. Bosons and fermions associated with a single factor 
U{k)mi,m2 of U{k;)^^ are interpreted as site variables; in an appropriate basis, the 
other surviving bosons and fermions are charged with respect to 2 factors and are 
therefore interpreted as link variables. 



3 Fermion action 



Here we examine the fermion determinant for the daughter theory. Expanded about 
the chosen point in moduli space, it is the fermion determinant for the lattice theory. 
We find it convenient to define 



Tr (Tf^T^TP) = Nr^P, 



''m.n 



"m.n'' 



(6) 



where A^ is an overall normalization that may be chosen as seems convenient, and 
m = (1711,1712) labels sites on a 2d square lattice; i = (1,0) and j = (0, 1) are unit 
vectors in the two directions. 

The fermion matrix depends on bosons 



^m ^m-' 



(7) 



as well as conjugates x^^ = [x^^Y , etc. The fermion action can be written in the form 
NV2 



(V^f,m , V^2^,m , €,m , " M^^^ 



52, n 

53, n 



The elements of the fermion matrix that follow from the expressions of [B] and the 
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conventions chosen here are: 



(M, 

(M, 
(M, 

(MZn] 

(Mz^: 



1,1 



11,2 



/ip ^ 
MP ' 



1,4 



2,3 



;3,i 



3,4 



4,2 



fP ■ ^ - 
m,n n+i 



^mn/2,2 



= (mZm 



.pup -V _ .^r- 
''m,n'^n m.n— x m) 



+pup V 
''m,n"''n+J 



''m,n+i-^n5 



ppv 

1- 

ppv 
m,r 

''m,nyn ^ ''m,n+ji'n5 
''m,n'^n ''m,n^n' 
''m,ny n 



mn)4,4 



= (M 

(MZ)i,3 

(MZ)2,1 
(M, 
TMP 



''m,nyn+i ''m,n+jy 



,pvp V 
m,n n+J 



''m.n^n? 



mn)2,4 



m,n— 7"m' 



m,n— 1 m' 



(M, 



mn)4,3 



3.1 Fermion zeromode 

The mother theory fermion modes A = A"T° that appear in Q are U{l)i ® U{1)2 
neutral; that is, they have ri = r2 = in (jSJ. As a consequence, the surviving parts 
correspond to T° that are nothing but the Cartan subalgebra of the mother theory, 
including the operator T° = 1. Because of the commutator in this fermion mode 
A° disappears from the action. Thus it corresponds to an ever-present zeromode in 
both the mother theory and the daughter theory. 

This zeromode eigenvalue of the daughter theory can be factored out following 
the method used in ^]. We deform the fermion matrix appearing in (jHI) according 
to 

M ^ M, = M + elTv^ (10) 

where Nj = Ak'^N'^ is the dimensionality of the fermion matrix and e <^ 1 is a 
deformation parameter that we will eventually take to zero. We factor out the zero 
mode through the definition 

M(0) = lim M(e), M(e) = e'^'^fM, det M(0) = lim e'Met . (11) 

If this deformation is added to the action, it explicitly breaks the exact lattice 
supersymmetry and gauge invariance. This infrared regulator could be removed in 
the continuum limit, say, by taking ea ^ A^~^. Noting that L = Na is the physical 
size of the lattice, the equivalent requirement is that e ^ L^^ be maintained as 
a — > 0, for fixed L. Thus in the thermodynamic limit {L oo), the deformation is 
removed. The parameter e is a soft infrared regulating mass, and is quite analogous 
to the soft mass introduced by Cohen et al. [H] in their Eq. (1.2) to control the 
bosonic zeromode of the theory.^ In the same way that the deformation introduced 
with fi does not modify the final results of the renormalization analysis of Section 

•^See also (HHJ below. 
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3.4 of [B], our e does not modify the result of the quantum continuum hmit. The 
essence of the argument is that we have introduced a vertex that will be proportional 
to the dimensionless quantity (72^a^ <^ g^a? /L, where g2 is the 2d coupling constant. 
Such contributions to the operator coefficients C, C in Eq. (3.29) of [H] vanish in 
the thermodynamic limit. Because the target theory is super-renormalizable, we are 
assured that the perturbative power counting arguments are reliable and the correct 
continuum limit is obtained. 

We have studied the convergence of detM(e) — >• detM(O). Indeed, we find that 
the convergence is rapid and that a reliable estimate for det M(0) can be obtained in 
this way. As a check, we have computed the eigenvalues of the undeformed matrix M, 
using the math package Maple, for a subset of 10 of the random boson configurations 
studied below. We find that the product of nonzero eigenvalues agrees with det M(0) 
in magnitude and phase to within at least 5 significant digits in each of the 10 cases. 

3.2 U{2) fermion determinant 

In our numerical work, we specialize to U{2). In notation introduced above, we choose 
T'^ = (I2, a"), Tr {T^'T'^TP) = 21^""" 

^000 ^ ^aOO ^ ^abO ^ ^ab ^ ^abc ^ -^abc ^^2) 

where underlining of indices indicates that all permutations are to be taken. 
The bosons have components 

Xm = xl^U + x^^a", = l/ml2 + Z/V> Zm = z'^U + z^a"" . (13) 

The lattice theory with lattice spacing a is obtained by expansion about a particular 
point in moduli space: 

av2 av2 

where ■ ■ ■ represent the quantum fiuctuations and all other bosons are expanded 
about the orgin. For this reason, in our study of det M(0) we scan over a Gaussian 
distribution where x'^,y^ have a nonzero mean l/a-\/2 = 1. The remainder of the 
bosons are drawn with mean zero. All bosons are taken from distributions with unit 
variance. 

For a set of 5000 draws on the bosons of a 2 x 2 lattice, we have extrapolated to 
e — i> and binned = argdet M(0) over its range, with 20 bins of size vr/lO. In Fig.^ 
we show the frequency for each bin, as a fraction of the total number of draws. It can 
be seen that once the zeromode eigenvalue is factored out, the product of the nonzero 
eigenvalues has arbitrary phase and that within statistical errors the distribution 
F{4>) is uniform: F(0) pa 1/20. Consequently, the effective lattice action [cf. ((T))] 
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Figure 1: Average frequency distribution for = argdet M(0), for 5000 random 
(Gaussian) draws, binned into intervals of vr/lO. Data was arranged into 50 blocks 
of 100 draws to estimate errors. The distribution of is seen to be, within errors, 
uniform. These results are for the U{2) lattice theory, with 2x2 lattice. 



is complex, with the phase of the fermion determinant a field-dependent quantity. 
This is to be contrasted with the observed reality of the fermion determinant in the 
mother theory [SB]- Presumably the orbifold projection does not commute with the 
conjugation operator that guarantees this reality in the mother theory. Indeed, it is 
not difficult to check that the orbifold projection operator does not commute with 
the usual conjugation operator that is involved in establishing the hermiticity of the 
A/" = 1 6d Minkowski spacetime action. Thus it is not surprising that the fermion 
determinant in the daughter theory is not real. 



4 Re-weighting 

It is worthwhile to explore whether or not the complex phase can be overcome for 
the purposes of simulation. A typical approach would be to compute averages of an 
operator O from the re-weighting identity: 

<^>^Wf '''' 

Here, = argdet M(0), as above, and "p.q." indicates phase- quenching: expectation 
values are computed with the replacement det M(0) — >■ | det M(0)|. Thus the effective 
bosonic action 

^p.q. = -ln|detM(0)| (16) 
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is used to generate the phase- quenched ensemble by standard Monte Carlo techniques. 
However, it is well-known that this tends to suffer efficiency problems: the number 
of configurations required to get an accurate estimate for, say, (exp(i0))pg, grows 
like exp(AF ■ Nj). Here, AF is the difference in free energy densities between the 
full ensemble and the phase-quenched ensemble. Recall that Nf is the dimensionality 
of the fermion matrix. It has been recently suggested how difficulties with this re- 
weighting approach might be surmounted by distribution factorization techniques 
Indeed, this method has been fruitfully applied in some other systems with complex 
action, some of which are quite similar to that studied here 123 12Z1 12H] • It would be 
interesting to see what progress may be made in the present context by applying this 
method. We are currently examining this idea and hope to report on it in a later 
publication. 

However, it is also possible that in the continuum limit the phase- quenched dis- 
tribution is sharply peaked at a value (po of (j). It follows that 



e 



i(f>0 



Thus it may be that the phase-quenched ensemble gives a good estimation in the 
continuum limit. For this reason it is of interest to study the distribution of (p as 
determined by the phase-quenched ensemble as a function of the lattice spacing. We 
next report results of a preliminary study of this distribution using Monte Carlo 
techniques. 



5 Phase-quenched distribution 

We generate a set of configurations, updating using the Metropolis algorithm applied 
to (|T?)|l . In this way, we sample the phase- quenched ensemble and estimate the cor- 
responding distribution of Fp q [(j)). For this purpose we need the bosonic action of 
the daughter theory. We specialize again to the U (2) case and introduce the notation 
Tr (Tt^T^T^T^) = 2t'^'^P^. Then refering to [HI, the bosonic action is given by 



S 



B 



^ n 



2 V-^n— i-^n— ? 



•^n-^n ~ i/n-Ji/n-j yni/n ~ ^n^n 



^n^n) 



-^n-^n ~^ Vn-jUn-j ynVn "I" ^n^n 



■P tjtA 



x'^x'fl + const. 
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with an implied sum over repeated superscripts. The parameter /i is a soft super- 
symmetry breaking mass that is inserted to stabihze the vacuum expectation values 
(jl4|) and is tuned to zero as the infinite volume limit is taken. In particular, in our 
simulations we are careful to respect the relation 

ga^ < fiaN < 1. (19) 

Further details on this deformation of the bosonic action may be found in jB]. 

We have studied the distribution Fp g {(f)) on small lattices for a few choices of the 
bare parameters. It is convenient to rewrite these parameters in terms of physically 
meaningful quantities with dimensions of length: 

a, g^' = {ga)-\ L = Na, (20) 

We expect that (7^^ is a rough measure of the correlation length for the system. Thus, 
we anticipate that systematics due to latticization and working at finite volume are 
kept to a minimum provided 

a < < L. (21) 

The effects of supersymmetry breaking in the infrared are expected to be negligi- 
ble provided g2^ <^ In the our simulations, we have examined the following 
parameter sets: 

{N;a,g^\L,fi-') = (2; 1/8, 8, 1/4, 1/2), (2; 1/20, 40, 1/10, 4/3), 

(3; 1/20, 40, 3/20, 2), (3; 2/3, 1, 2, 2), 
(4; 1/2, 1,2, 2). (22) 

The first three sets suffer from extreme finite volume effects, but the effects of dis- 
cretization are expected to be rather small since ag2 < 1/64. Supersymmetry breaking 
is expected to be small in the second and third sets since fi/g2 < 3/160. The last two 
cases respect a weaker version of the "physical" constraint (EH), namely a < (7^^ < L, 
and simultaneously a weaker version of (fT^ . namely ga"^ < fiaN = 1. (The condi- 
tions were of necessity weakened, due to the small lattices we are presently working 
with.) In all cases, we find that the distribution is uniform; that is, up to statistical 
fiuctuations, the results are indistinguishable from those of Fig. [T] 

6 Conclusions 

Naturally we would like to study the phase- quenched distribution Fp.q.(0) on larger 
lattices where we can better approach the continuum limit. However, these simula- 
tions are rather demanding since they involve a fermion determinant; we leave such 
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an effort to future work, as it will take some time to accumulate the necessary lattice 
data and develop efficient algorithms. Moreover, it is not at all clear that the phase of 
the determinant over all lattice fermion modes, over all nonzero measure boson con- 
figurations in the phase-quenched ensemble, is entirely relevant to the phase of the 
continuum limit fermion determinant. Only those modes that have a significant over- 
lap with the low energy interpolating fields are expected to have physical meaning. 
The other modes integrate out without significantly affecting correlation functions of 
physical operators. This is closely related to the distribution factorization method 
mentioned above, in the case where the operator O appearing in ()15|) corresponds to 
a good approximation to an operator in the target theory with external momentum 
scales well below a~^. Research in this direction is underway. 

Alternatively, it may be hoped that a different orbifold may preserve the reflection 
positivity of the mother theory while still producing an effective lattice theory with 
exact supersymmetry. We suspect that the fermion determinant would be positive in 
such a construction. In the present constuction it seems that reflection positivity is 
only recovered in the continuum limit. We are also pursuing this notion in ongoing 
research. 

The complex action problem encountered here provides an amusing opportunity. 
Some exact results are available for the (4,4) 2d super- Yang-Mills theories. In partic- 
ular, the 1-loop effective action has been determined j29j, modulo higher derivative 
terms. It would be interesting to compare the results of complex action simulation 
techniques to analytic results such as these. ^ It may be that the supersymmetric 
lattice actions proposed here provide an independent check on Monte Carlo methods 
that have been suggested in the simulation of QCD at finite temperature and baryon 
density. 
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